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Abstract 

We present a consistent picture of tunnelling in field theory. Our results apply both 
to high-temperature field theories in four dimensions and to zero-temperature three- 
dimensional ones. Our approach is based on the notion of a coarse-grained potential 
Uk that incorporates the effect of fluctuations with characteristic momenta above a 
given scale k. Uk is non-convex and becomes equal to the convex effective potential 
for A; ^ 0. We demonstrate that a consistent calculation of the nucleation rate must 
be performed at non-zero values of k, larger than the typical scale of the saddle- 
point configuration that dominates tunnelling. The nucleation rate is exponentially 
suppressed by the action Sk of this saddle point. The pre-exponential factor Ak, 
which includes the fluctuation determinant around the saddle-point configuration, 
is well-defined and finite. Both Sk and Ak are /c-dependent, but this dependence 
cancels in the expression for the nucleation rate. This picture breaks down in the 
limit of very weakly first-order phase transitions, for which the pre-exponential factor 
compensates the exponential suppression. 



1 Introduction 



The consistent description of first-order phase transitions is a difficult problem which has attracted the 
attention of statistical and particle physicists for a long time (for a review see ref. and references 
therein) . Our present understanding of these phenomena is based on the work of Langer on nucleation 
theory [||. His formalism has been applied to relativistic field theory by Coleman Q and Callan Q 
and extended by Affleck and Linde |^] to finite temperature. The basic quantity in this approach is 
the nucleation rate, which gives the probability per unit time and volume to nucleate a certain region 
of the stable phase (the true vacuum) within the metastable phase (the false vacuum). The rate 
is exponentially suppressed by the free energy of the critical bubble, which is a static configuration 
(usually assumed to be spherically symmetric) within the metastable phase whose interior consists of 
the stable phase. This configuration has a certain radius that can be determined from the parameters 
of the underlying theory. Bubbles slightly larger than the critical one expand rapidly, thus converting 
the metastable phase into the stable one. Possible deformations of the critical bubble generate a static 
pre-exponential factor in the nucleation rate. The leading contribution to this factor has the form of a 
fluctuation determinant. Another dynamical prefactor determines the fast growth rate of the bubbles 
that are slightly larger than the critical one (2|, ^ . In this work we concentrate on the calculation of the 
static prefactor. We are, therefore, mostly concerned with the rate of nucleation of critical bubbles. 
Their real-time evolution after their nucleation is a separate question, which we hope to address in a 
future publication. 

The nucleation rate per unit volume / (probability of nucleation of a critical bubble per unit 
time and volume) for a four-dimensional field theory at temperature T, in the limit that thermal 
fluctuations dominate over quantum fluctuations, is given by 



2tt \2tt) 



exp(-5). (l.I) 



Here T is the free energy of the system for a given configuration of the field </> that acts as the 
order parameter of the problem. The rescaled free energy of the critical bubble is S = T^/T = 
[r((^b('^)) — r(0)] /T, where (ph{r) is the spherically-symmetric bubble configuration and (f) = cor- 
responds to the false vacuum. The fluctuation determinants are evaluated either at cf) = or around 
(f> = (phif)- The prime in the fluctuation determinant around the bubble denotes that the three zero 
eigenvalues of the operator [S'^F / 6(p'^](f,=^^ have been removed. Their contribution generates the fac- 
tor (5/27r)'^/^ and the volume factor that is absorbed in the definition of / (nucleation rate per unit 
volume). The quantity Eq is the square root of the absolute value of the unique negative eigenvalue. 
This last contribution appears only for the high-temperature theory |^]. It is absent in the expression 
for the quantum-tunnelling rate in the zero-temperature three-dimensional theory. 

The free energy of the critical bubble can be easily determined either analytically or numerically. 
The bubble is the dominant saddle-point configuration that interpolates between the two vacua. Its 
profile is determined by a differential equation, which can be integrated numerically or even solved 
analytically in simple cases. The calculation of the fluctuation determinants in the prefactor is a more 
difficult task. They can be brought in a more manageable form if one employs spherical coordinates. 
However, the true difficulty concerns the ultraviolet divergences that are inherent in their calculation. 
An appropriate regularization scheme must be employed in order to control them |§-|1C]. 



The situation becomes even more complicated in the case of radiatively induced first-order phase 
transitions. These are a consequence of the appearance of a new vacuum state in the theory as a 



result of the integration of (quantum or thermal) fluctuations |11]. In field theory the free energy 
(more precisely the thermodynamic potential) density of a system for homogeneous configurations is 
usually identified with the temperature-dependent effective potential. A radiatively induced first-order 
phase transition appears in theories for which the tree-level potential has only one minimum, while a 
second minimum appears at the level of radiative corrections, usually computed within a perturbative 
scheme. This approach, however, faces two fundamental difficulties: 
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a) The effective potential, being the Legendre transform of the generating functional for the con- 
nected Green functions, is a convex function of the field. Consequently, it does not seem to be 
the appropriate quantity for the study of tunnelling, as no structure with more than one minima 
separated by a barrier existsQ. 

b) The fluctuation determinants in the expression for the nucleation rate have a form completely 
analogous to the one-loop correction to the potential. The question of double-counting the effect 
of fluctuations (in the potential and the prefactor) must be properly addressed. 

In this paper we demonstrate that all the above issues can be resolved through the implemention of the 
notion of coarse graining in the formalism. The appropriate quantity for the description of the physical 
system is the effective average action , which is the generalization in the continuum of the blockspin 
action of Kadanoff H^. The dependence of this action on the coarse-graining scale k is described by 
an exact flow equation [l^ typical of the Wilson approach to the renormalization group The 
formalism has been applied with success to second-order phase transitions. A complete picture has 
emerged for the phase transitions in a variety of scalar models, with a reliable determination of both 
non-universal quantities (such as critical temperatures) and universal ones (such as critical exponents, 
the equation of state and crossover curves) [|r7|-p^. The generalization of the formalism to gauge 



theories |23| has led to the study of second-order and first-order phase transitions for the Abelian and 



non-Abelian Higgs models |^^-[^], with implications for the electroweak phase transition |27]. The 



method has also been applied to the chiral phase transition |28|. The framework for the discussion of 



nucleation rates in first-order phase transitions has been set in ref.s |^, 3C]^. 

We first summarize the basic notions in the calculation of nucleation rates in our approach. We 
define the bare theory at some high scale A that can be identified with the ultraviolet cutoff. The 
renormalized theory at lower scales k is described in terms of the effective average action F^, which can 
be interpreted as the coarse-grained free energy at a given scale k. Fluctuations with characteristic 
momenta ^ /c^ are integrated out and their effect is incorporated in r^. The k dependence of 
Ffc is determined by an exact flow equation. This can be translated into evolution equations for the 
invariants appearing in a derivative expansion of the action. We consider only the effective average 
potential and a standard kinetic term and neglect higher derivative terms in the action. The 
validity of our assumption is discussed in the next section. At scales k below the temperature T, the 



theory can be described in terms of an effective three-dimensional action at zero temperature |17, 27|. 
This dimensional reduction indicates the absence of explicit time dependence for the parameters of the 
theory. It is a consequence of our implicit assumption that the high-frequency modes of the system 
are in thermal equilibrium and their time dependence can be averaged out. On the contrary, the 
low-frequency modes have real-time dynamics that are related to the behaviour of the system during 
and after nucleation. In this work we concentrate on static properties of the system, such as the 
characteristics of the critical bubble and the nucleation rate. For such quantities, the description in 
terms of an effective three-dimensional theory is sufficient. The study of dynamical questions, such 
as the expansion rate of bubbles slightly larger than the critical one, requires a study of the full 
four-dimensional theory, at least for the low- frequency modes. 

We determine the form of the potential Uk at scales k below the temperature through a numerical 
solution of the evolution equation. We are interested in theories for which has two minima separated 
by a barrier for low values of k. This structure may exist already at the level of the bare potential, 
or appear as the result of the integration of the high-frequency modes. When k"^ becomes smaller 



^ It has been argued in ref. that the appropriate quantity for the study of tunnelhng is the generating functional 
of the IPI Green functions (calculated perturbatively), which differs from the effective potential in the non-convex 
regions. However, as we show in the following, the consistent picture must rely on the notion of coarse graining and 
on the separation of the high-frequency fluctuations that are responsible for the non-convexity of the potential, from 
the low-frequency ones that are relevant for tunnelling. Such notions cannot be easily implemented in the context of 
perturbation theory. 

^For a related work see ref. |3l|. 
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than the typical positive curvature of the convex parts of C/^, the massive modes that induce the 
evolution of the potential decouple. As a result the convex regions of the potential stop evolving. 
However, the non-convex part (the barrier) continues to evolve. Full convexity is approached in the 



limit /c — > |l^, 2C]. It is a consequence of the integration of configurations that interpolate 
between the minima in the functional integral that defines [p2| , |3^ . All explicit information about 
tunnelling is lost in the resulting Maxwell construction for the effective potential. 

It is clear from the above that the calculation of the nucleation rate must be performed at a 
non-zero value of k, such that configurations that interpolate between the minima are not taken into 
account. This value must be chosen so that the convex parts of the potential have stopped evolving 
significantly, while a well-defined barrier still exists. For a range of values of k that satisfy this 
requirement we perform the calculation of the nucleation rate. The profile and the free energy of the 
bubble are determined in the standard way. The evaluation of the fiuctuation determinants in the 
prefactor is again performed following standard techniques, but with an important modification. An 
ultraviolet cutoff of order k is imposed, such that fluctuations with characteristic momenta q'^ ^ k"^ 
are not included. The reason is that the effect of such fluctuations is already incorporated in Uk- This 
modiflcation resolves two of the serious problems mentioned earlier. The pre-exponential factor is now 
finite and no double-counting of the fluctuations takes place. 

However, an important issue arises at this point. The scale k was introduced in the problem as 
a mere calculational tool. If our approach makes sense, the choice of k should not affect physical 
parameters such as the nucleation rate. The remarkable outcome of our study is that this expectation 
is conflrmed. Despite the significant k dependence of the free energy of the bubble and the prefactor, 
the nucleation rate is fc-independent to a good accuracy. A residual small k dependence can be 
interpreted as a measure of the contribution of the next order in the saddle-point approximation for 
fluctuations around the bubble. 

In the following sections we present the details of the calculation outlined above for a theory of a 
real scalar fleld. In order not to obscure the essential physics by complicating the model too much, 
we have made some simplifications. We do not discuss the evolution of F^ for k T. For readers 
who are interested in the details of the mechanism of dimensional reduction in our approach, detailed 



discussions can be found in ref.s |17, 22, 27| for a variety of models. We start the evolution at a scale 
ko sufficiently below the temperature of the system, so that the dynamics is three-dimensional to a 
good approximation. As an initial condition we consider a potential U^q with two inequivalent minima 
separated by a barrier. This form of the potential is determined by the bare potential U/\_ and the 
integration of fluctuations between the scales A and k^. Some of these fluctuations may correspond 
to additional massive degrees of freedom that decoupled above the scale k^. In the next section, we 
describe in detail the initial form of the potential we use. We integrate the evolution equation for 
the effective three-dimensional theory starting at the scale k^, and perform the calculation of the 
nucleation rate as described above. A signiflcant evolution, with a substantial variation of the form 
of the potential, may take place between ko and the range of scales where we compute the nucleation 
rate. 

In the following section we derive the evolution equation for the potential and describe the initial 
condtion for its integration. In section 3 we discuss the formalism we employ for the calculation of 
the nucleation rate. Our results are presented in section 4 and our conclusions are given in section 5. 



2 Evolution equation for the potential 

In this section we summarize the formalism of the effective average action for a theory of a real scalar 
field (p and derive the evolution equation for the potential. We discuss the effective three-dimensional 
theory that results from the dimensional reduction of a high-temperature four-dimensional theory at 
scales below the temperature. The temperature can be absorbed in a redefinition of the field and its 
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potential, so that these have dimensions appropriate for an effective three-dimensional theory 



U{^) = (2.1) 

In this way, the temperature does not appear explicitly in our expressions. This has the additional 
advantage of permitting the straightforward application of our results to the problem of quantum 
tunnelling in a three-dimensional theory at zero temperature. 

2.1 Intuitive derivation 

Before presenting the rigorous derivation, it is instructive to derive the evolution equation for the 
potential based on an intuitive argument, along the lines of ref.s |13, We start by considering the 



Z2-symmetric scalar model, in Euclidean three-dimensional space. The one-loop effective potential is 
given by the expression 



U^\p) = V{p) + hndet[Pk + V'{p) + 2V"{p)p\ 



= V{p) + i 1^ ^ In [P,iq) + V'ip) + 2V"ip)p] , (2.2) 

where V{p) is the bare potential. In order to be consistent with the conventions in previous publica- 
tions, we have defined the variable 

P = \<t^\ (2.3) 

which we frequently use in this section. Primes denote derivatives with respect to p: V'{p) = dV /dp. 
In terms of this variable, the mass term of the scalar field is d'^V/d(j)^ = V'{p) + 2V"{p)p. The inverse 
propagator Pk{q) in momentum space for a massless field is given by Pk=o{q) = in perturbation 
theory. We assume that the momentum integration is regulated by an ultraviolet cutoff A. 

We would like to introduce an effective infrared cutoff k for the low-frequency modes, so that 
the momentum integration in eq. (^]^) does not receive contributions from modes with characteristic 
momenta q'^ ^ k'^. The simplest way to achieve this is through the addition of a mass term k'^ to the 
perturbative inverse propagator, so that for a massless field 

Pk{q)=q^ + k\ (2.4) 

The potential now depends on k, as indicated by the subscript in eq. (p.2|). 

The next step is to derive an evolution equation for the change of C/^ with the scale k and follow 
the evolution for k ^ 0. For this purpose we take the logarithmic derivative with respect to k and 
substitute Uk for V in the right-hand side of eq. (^.2|). The intuitive justification for this replacement 
is based on the fact that the new contributions to the momentum integration, when k is lowered by 
a small amount Ak, come from the region k — Ak < q < k. The relevant mass term and couplings 
that should appear in the evolution equation are the renormalized ones at the scale k (which, for 
the scalar field, are related to derivatives of Uk) and not the bare ones. This "renormalization-group 
improvement" results in the evolution equation 

dU'kip) ^ 1 f d^q dPk W'^{p)+2U'^{p)p 
9t 2 7 (27r)3 dt [p,(g) + c/^(^) + 2^^(p)p]2 

e 3U'^ip) + 2U'^'{p)p 

8^ ./k^ + U',ip) + 2U'^{p)p 
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where t = ln(fc/A), with A identified with the ultraviolet cutoff of the theory. We have derived the 
evolution equation for Ul^{p) (and not Uk{p)) because this is the easiest to integrate numerically. For 
k = A the infrared and ultraviolet cutoffs coincide, and no integration of fluctuations takes place. This 
determines the initial conditions for the solution of eq. ( |2.6D as Ux{p) = ^(p)- In the opposite limit, 
/c — > 0, one recovers the effective potential U{p) = Uo{p). 

One could look for an iterative solution of eq. ( p.5| ) . The first iteration results in the equation 

2 det [Pk + U^{p) + 2U^ {p)p\ 



for the effective potential. This expression can be compared with eq. ( |2.2[ ) with k = 0. The two 
equations have the same structure, but the bare potential V{p) is replaced by the fe-dependent po- 
tential Ukip)- Also, the determinant resulting from the radiative corrections is replaced by a ratio 



of determinants. The numerator of this ratio is what one would expect from eq. (2^) with k = 0. 
The denominator involves the inverse propagator P^. and effectively removes the fluctuations with 
characteristic momenta <^ k'^. This is justified by the fact that the effect of these fluctuations has 
already been incorporated in Uk{p)- There is no need for the addition of an extra regulator A as in 
eq. ( |2.2| ). Its role is played by k, which acts as an effective ultraviolet cutoff in the calculation of the 
effective potential from Uk{p)- 

The above observation will be important in the following section where nucleation rates will be 
computed. The first part of our procedure will be to integrate the evolution equation ( |2.5D numeri- 
cally, from a scale /cq, where we choose the initial form of the potential Uk^ip)^ down to a non-zero 
scale k. This effective integration of the high-frequency modes generates the appropriate potential 
that describes the dynamics of low- frequency modes with ^ k'^. Tunnelling will be discussed in 
the context of the low-energy theory through the standard saddle-point approximation [^, ^. The 
first correction to the leading semiclassical result involves the usual "one-loop" form of a fluctuation 
determinant det [Pk=o + O], where O is an operator to be defined in the next section. A consistent 
treatment of the effect of fluctuations, which avoids double-counting the high-frequency modes, can 
be obtained if this correction is modified to det [{Pk=o + O) / {Pk + O)] as in eq. (^). 

2.2 Rigorous derivation 

The evolution equation for the potential can be derived within a more rigorous approach through 



the formalism of the effective average action |13, 15, The effective average action F^, for a 



theory described by a bare action 5, results from the effective integration of degrees of freedom 
with characteristic momenta larger than a given infrared cutoff k. Its dependence on the scale k is 
described by an exact flow equation. In this subsection we summarize the formalism for the case of a 
Z2-symmetric theory of a real scalar field in Euclidean three-dimensional space. A detailed discussion 
can be found in ref.s jl^, |l8|. 

We specify the action together with some ultraviolet cutoff A, so that the theory is properly 
regulated. We add to the kinetic term a piece that has the following form in momentum space 



^kS[x] = ljd\ Rk{q)x*{q)xiq), (2. 



where x*{q) = xi~Q)- The function Rk is used to prevent the propagation of modes with characteristic 
momenta ^ k"^. As a result, the inverse propagator for the action S + AS has a minimum ~ fc^. 
The modes with (f' ^> k"^ are unaffected, while the low-frequency modes with q^ <C k"^ are cut off: 

lim Rk^k^. (2.9) 
We emphasize at this point that many alternative choices of Rk are possible. 
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We subsequently introduce sources and define the generating functional for the connected Green 
functions for the action S + AS. Through a Legendre transformation we obtain the generating func- 
tional for the IPX Green functions [</>] , where <j) is the expectation value of the field x i^i the presence 
of sources. The use of the modified propagator for the calculation of Tk results in the effective in- 
tegration of only the fluctuations with q'^ ^ k'^. Finally, the effective average action is obtained by 
removing the infrared cutoff 

rk[4>] = f fe[0] - ^ 1 d\Rk{qW{q)m. (2.10) 

For k equal to the ultraviolet cutoff A ^, F^ becomes equal to the bare action S (no effective 
integration of modes takes place) , while for A; — > it tends towards the effective action F corresponding 
to S (all the modes are included). The interpolation of Tk between the bare and the effective action 
makes it a very useful field-theoretical tool. The means for practical calculations is provided by an 
exact flow equation]^, which describes the response of the effective average action to variations of the 
infrared cutoff {t = ln(A;/A)) [|l|]: 

^T,m = -Tr {(Ff^M + i?,)-i|-i?fcl . (2.11) 



dt 2 L-rj ^; 

(2) 

Here F^ is the second functional derivative of F^ with respect to 0. 

Making use of the Z2 symmetry, we parametrize the effective average action as 

Tk = j d^x S^Ukip) + ^5^'/' -Zkip, -d^y. d^cP (2.12) 

where the normal ordering indicates that the derivative operators are always on the right. The dots 
stand for invariants that involve more derivatives of the field. In order to turn the flow equation for 
the effective average action into equations for Uk, Z^, etc, we have to evaluate the trace in eq. ( |2.11| ) 
for properly chosen background field configurations. For the evolution equation for Uk we have to 
expand around a constant field configuration. We find [l^, |l8| 

l/J^^ dRkiq)/dt 

dt'^P' 2J (2^)3 Zk{p,q^)q^ + Rk{q) + U',{p) + 2pU'f!{py ^ ' 

The above equation is an exact evolution equation for the potential. Its solution, however, re- 
quires information on the wave- function renormalization Zk{p,q'^)- Throughout this paper we shall 



set Zk{p,q^) = 1, which corresponds to the first order of the derivative expansion of eq. ( 2.12 ). This 
is expected to be a good approximation because the size of Zk is related to the anomalous dimension 
of the field, which is small for this model (r/ ~ 0.035 for the three-dimensional theory we consider). 
For r/ = the kinetic term in the /c-dependent inverse propagator must be exactly proportional to q^ 
both for — > and q^ ^ co. Several studies have confirmed the smallness of the corrections arising 
from the deviation of Zk{p,q'^) from 1 |18, 35, 21|. Within our approximation, eq. ( |2.13| ) reproduces 
eq. (|2.5|) through the definition 



Pk{q)=q' + Rk{q). (2.14) 

For the choice 

Rk{q) = k\ (2.15) 

eq. (|2.6| ) is obtained. 

^ For scales k > T the high-temperature behaviour of the four- dimensional theory is relevant. For this reason our 
discussion in terms of an effective three-dimensional theory is not sufficient. In this work, however, we never discuss the 
evolution at such high scales, but work at scales k < ko ^ T instead. We refer to the ultraviolet cutoff A only at this 
point for reasons of completeness of the presentation. For a full discussion of the evolution at large momentum scales 
and the mechanism of dimensional reduction, see ref.s p7| . 

* For other versions of exact renormalization group equations, see ref.s jl^, 
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We should point out that the choice of eq. ( p. 15 ) for the cutoff function Rk{q) may be problematic. 



For example, in four dimensions it results in divergent integrals in the right-hand side of evolution 
equations such as eq. ( |2.13| ). Even for the three-dimensional case that we are considering, the right- 
hand side of eq. ( |2.13| ) involves an irrelevant p-independent divergent constant, which disappears in 
the evolution equation for f7^(/o). For these reasons, in most cases it is preferable to work with a cutoff 
function of the form 

1 — exp (— g^/K^j 

for which integrals such as the one in the right-hand side of eq. ( |2.13| ) are finite. However, such a cutoff 
function corresponds to a non-local operator in position space, which would make the calculations of 
the next section impossible. Because of this, we employ the cutoff function of eq. ( |2.15| ) in this 
paper and emphasize that care must be taken in the extension of our discussion to the case of full 
four-dimensional dynamics. 

2.3 Scale-invariant form of the evolution equation 

It is convenient to cast the evolution equation (2.6) in a form that does not explicitly depend on the 



scale k. This makes the identification of possible fixed points easier. For this reason we define the 
dimensionless quantities 



Uk{p) 



Ukjp) 



p = (2.17) 

Primes on uj. denote derivatives with respect to p. We can now rewrite the evolution equation for the 
potential as 

du' 1 

_A= _24 + p4^ + _(3< + 2K')i?K + 2<p). (2.18) 

The non-trivial solution of the above equation with du'^/dt = corresponds to the Wilson-Fisher 
fixed point that determines the dynamics of the second-order phase transition in the Z2-symmetric 



theory |l|, |21 ] . 



The dimensionless function L\{w) is given by 

1 f dPk 



= --/ dxV^^{Pk + w)-\ (2.19) 
k Jo ot 

with X = . This has been discussed extensively in ref.s |13, 18, |3^ (for various forms of the infrared- 
regulating function Rk and in various dimensions). It has the interesting property that it falls off 
for large values of w following a power law. As a result it introduces a threshold behaviour for the 
contributions of massive modes to the evolution equation. The third term in the right-hand side of 
the evolution equation ( |2.18| ) includes the L\ function with the mass of the cj) field divided by k"^ as 
its argument. When the scale k"^ crosses below the running squared mass this contribution vanishes 
and the massive mode decouples. As a result the evolution of Uk{p) stops. The function L\{w) also 
has a pole at w = —1. This property induces the convexity of the potential in the limit — > 0. The 
argument of Li{w) in the non-convex regions is given by the negative curvature of the potential divided 
by k'^. As the pole cannot be crossed, the curvature follows k to zero, thus inducing the convexity of 
the effective potential. For the choice of eq. ( p. 15 ) for the cutoff function, we obtain 



Lliw) = -^=^, (2.20) 
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in agreement with eq. ( p7 

Two algorithms for the numerical integration of eq. ( |2.18| ) have been presented in detail in ref. [35|. 



The comparison of the two methods provides a good check on possible systematic numerical errors. 
The two algorithms give results that agree at the 0.3% level. We expect the numerical solution to 
be an approximation of the solution of the partial differential equation (|2.18| ) with the same level of 
accuracy. 

2.4 Explicit breaking of the Z2 symmetry 

Up to this point we have not discussed the breaking of the Z2 symmetry that could lead to vacuum 
instability. The formalism of the previous subsections is completely Z2-invariant and has to be modified 
in order to account for the symmetry-breaking effects. However, we shall work within a framework for 
which the modifications are minimal. As we explained in the introduction, we consider an effective 
three-dimensional theory that results from the dimensional reduction of a high-temperature four- 
dimensional one. We define this theory at a scale ^ T. The form Uk is determined by the 
bare potential U\ and the integration of fluctuations between the scales A and /cq- Some of these 
fluctuations may correspond to additional massive degrees of freedom that decoupled above the scale 
ko- We choose a form of the potential that breaks the Z2 symmetry, while permitting the presence of 
a non-trivial evolution at scales k < ko. 

Following ref. |^6| , we consider theories that are described by potentials of the form 

Ukoi^) = ^mlj^ + ^jko^^ + ^Xkj\ (2.21) 



By a variable shift 



^ = <P + ^ (2.22) 



we can bring this potential to the form 

^fco(^) = Cfco - J^o- + ^^io^^ + ^Afcgcr^, (2.23a) 

with 

3 

T = .2.23b) 

6Afco 

The right-hand side of the exact flow equation for the effective average action of the a field is not 
affected by the linear term —J^a or the constant c^q. Therefore, the evolution equation for the 
potential can be integrated using the Z2-symmetric formalism. The potential of (p is recovered through 
the relation 

UkicP) = Ck, - J,a + U^\a) = c^, - • + ^) + f/f [c^ + ^) , (2.24) 
where U^'^ is the potential of the Z2-symmetric model with mass term ji^^ and quartic coupling A^g. 



In ref. [36| it has been verified that the above procedure gives the same result as the straightforward 



integration of the evolution equation ( p.l8| ) with an initial condition given by eq. (2.21). 



In the following we integrate numerically eq. ( |2.18| ) for Z2-symmetric models. Different types of 
evolution can be obtained by keeping A^^ constant and varying the mass term More specifically, 
the Wilson-Fisher fixed point of the three-dimensional theory can be approached for a certain (neg- 
ative) value (Aifc(j)cr- This fixed point determines the properties of the second-order phase transition 
of the Z2-symmetric theory. Through the appropriate choice of 7^;,, we can vary the difference in 
the energy density between the minima of the potential according to eq. ( |2.24| ). Arbitrarily weakly 
first-order phase transitions can be studied in the limit — >■ {^i\^)cr-, Iko 0. 
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3 Calculation of the nucleation rate 



3.1 General formalism 



As we discussed in the previous sections, we study tunnelling in a theory of a real scalar field with a 
Euclidean action 

Tk = I d^x + C/fe((/))| (3.1) 

and a coarse-graining scale k 0. The parameters in the above action are effective three-dimensional 
ones defined according to eqs. ( |2.1D . The potential Uk{(j)) has two minima corresponding to vacua with 
different vacuum energy densities: the stable (true) minimum is located at (j) = (j)t and the unstable 
(false) one at cp = cpj = 0. 

The unrenormalized decay rate per unit volume from the false minimum towards the true one is 
given by |, | | 



2tt Uvrj 



3/2 



det'[52r^ 



det [52^^/502 



l</>=0 



-1/2 



exp(-S'fc 



(3.2) 



This is analogous to eq. ( ]1.1| ) after the absorption of the explicit factors of T in the redefinition of the 
field and potential and the introduction of a coarse-graining scale. Notice that if the factor E^jlTx 
is removed, the above expression reproduces the quantum-tunnelling rate of the zero-temperature 
three-dimensional theory |Q, As this factor gives only a small contribution to the total rate, our 
discussion applies to quantum tunnelling as well within a good approximation. 

The nucleation rate is exponentially suppressed by the action S'fc (the rescaled free energy) of the 
bubble configuration ^\:,{r'). This is an S0(3)-invariant solution of the classical equations of motion 
which interpolates between the local maxima of the potential —Uk{4')- It satisfies the equation 



b _^ 2 d(j)^ 

dr"^ r dr 



(3.3) 



with the boundary conditions </>b ^ for r ^ oo and d(j)]^/dr = for r = 0. The action Sk of the 
bubble is given by 



rc 

Jo 



1 /#b(r) 



dr 



+ UkiMr)) 



r''dr^Sl + Sl, 



(3.4) 



where the kinetic and potential contributions, 5^ and respectively, satisfy 



SI 

Si 



1 

3" 



(3.5) 



The pre-exponential factor corresponds to the first correction to the semiclassical approximation 
in the saddle-point method. The numerator is the fluctuation determinant around the bubble 



det' 



det' \-d^ + Uj:{(t> = Mr)) 



(3.6) 



while the denominator is the fluctuation determinant around the false vacuum 



0; = O 



det 



0=0 



det 



(3.7) 



In this section we revert to the standard way of denoting the derivatives with respect to (p with primes. 
This should not be confusing as we always indicate the quantity with respect to which we differentiate 
as the argument of the function. The differential operator —d^ + (</'b(^)) has three zero modes (the 
three spatial translations of the bounce). The prime over the determinant indicates that these modes 
have to be omitted in its calculation. Their contribution generates the factor (5'fc/27r)^/2 gq (^) 
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and the volume factor that is absorbed in the definition of / (nucleation rate per unit volume). The 
quantity Eq is the square root of the absolute value of the unique negative eigenvalue of the above 
operator. 

The pre-exponential factor defined in eq. is in general ultraviolet-divergent and an appropriate 
regularization scheme must be employed. Within our approach, the form of the regularization is 
dictated by the discussion at the end of subsection 2.1. The effect of the high-frequency modes has 
been incorporated in the form of the coarse-grained potential Uk{4)), which is obtained through the 
integration of the evolution equation ( 2.1§| ). In order not to double-count this effect, fluctuation 
determinants computed within the low-energy theory must be replaced by a ratio of determinants, 
in complete analogy to eq. ( |2.7| ). This implies that a consistent expression for the nucleation rate is 
given by 



/ 



Akexp{-Sk) 

271 \2ttJ 

det' [-d^ + U'^{4,^{r))] 



det [-92 + (-^2) + [/^(O)] 



det [-52 + ^7^(0)] 



-1/2 



(3. 



The infrared-regulating function has been discussed in the previous section (see eq. ( 2.14| )). The 
simplest choice for this function, which is necessary for the feasibility of the computation of the next 
section, is given by eq. ( p. 15 ), i.e. Rk = k"^- We point out that only the operator —d^ + U'f!{(t)h{r)) has 
negative and zero eigenvalues that require special treatment in eq. (|3.8|). 



3.2 Details of the numerical computation 

Prom this point on, we use the form R^ = for the infrared-regulating function. The profile of the 
bubble can be easily computed with the "shooting" method [^. We integrate eq. (|3.3| ) numerically, 
starting at r = with a value of 4> near the true minimum (pt and d<j)/dr = 0. We then adjust the 
initial value of (p so that the boundary condition 4>i, ^ for r — > 00 is satisfied. (In practice this 
condition is satisfied at a value Too sufficiently larger than the typical size of the bubble.) 

The computation of the fluctuation determinants is more complicated. The differential operators 
that appear in eq. (|3.8[) have the general form 



where 



= -d^ + ml + aWkir) 

ml ^ U'a^) + Kk\ 
Wk{r) ^ U'^{Mr))-U'f!{0), 



(3.9a) 

(3.9b) 
(3.9c) 



with K, a = or 1. Since the W^a operators are S0(3) symmetric, it is convenient to use spherical 
coordinates {r,9,(p) and express the eigenfunctions ip in terms of spherical harmonics: 'ijj{r,6,(p) = 
u{r)/r Y£jn{0, (p) |§, |9|. Here i and m are the usual angular quantum numbers. The Laplacian operator 
(?2 takes the form 



52 



so that 



(i2 1(1 + 1) 



dr 



--v,V, 



(3.10) 



\2l+l 



detW«„ = Y[{detWe^^^, 

Wfca = -Vj + ml + aWk{r). 



(3.11) 
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Figure 1: The steps in the computation of the nucleation rate for a model with li^^Jk^ = —0.05, 
•^feo/^o = 0-1; Iko/ky^ = —0.0634. The dimensionful quantities are given in units of kf = 0.223 /cq. 



We recall that det Wi^a is defined as the product of all eigenvalues A that lead to solutions of 
WiKa'^i^) = ^u{r), with the function u(r) vanishing at r = and r oo. The computation of 
such complicated determinants is made possible by a powerful theorem ^ that relates ratios of 
determinants to solutions of ordinary differential equations. In particular, we have 



detWtei _ det[-V| + ml + l- Wkir)] _ ye^i{r 



oo 



det Wfco det[- V2 + ml + 0- Wk{r)] 
where yeK,a{f) is the solution of the differential equation 



yiKoir oo) ' 



(3.12) 



dr'^ 



+ 



+ ml + aWk{r) 



yiKair) = 0, 



(3.13) 



with the behaviour yi^ai''') f^^^ for ^ 0. Such equations can be easily solved numerically with 
Mathematica |39]. In the "free case" (a = 0) it is possible to obtain the exact analytical solution^ 



VIkQ oc ii+i/2{mnr)/y/m~r. 



(3.14) 



^The function i„ is the standard Bessel I function, defined as i„(z) = Bessell[n, z\ in Mathematica notation For 
n = £+1/2 a semi-integer, i„ can be expressed in terms of elementary functions as y/zie+-i/2iz) = P^{1/ z)e^ +P!l{1/ z)e^' , 
where P± are polynomials of degree £. For z oo (or, more precisely, for z ^ £) ig^xjiiz) — > e'^ /^2-kz. 
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The final expression for the nucleation rate, appropriate for an efficient numerical computation, is 



3/2 



Co 



2tt \2i: J 
V Ifool 



exp {-Sk) cqCi q, 

£=2 



1/2 



Cl 



9m) 



(2£+l)/2 



(3.15) 



The factors q for l>2 can be computed in a straightforward way through eqs. ( 3.12| ), ( |3.13| ) as we 
explained above. The calculation of ci is more complicated because of the necessity to eliminate the 
zero eigenvalues in g'^^. This can be achieved by replacing the operator Wioo with Wioo + e U'j^{0) 
and evaluating g'lQ as 



a'lo 



The (unique) negative eigenvalue —Eq of Wool can be obtained by solving the equation Wooi^* = —EqU 
and using the shooting method to determine the value of Eq that ensures the correct boundary 
condition u{r ^ oo) = 0. 

As a final remark, we give the explicit expression for q in the limit of large i. It can be obtained 
by solving the differential equations ( 3.13| ), using first-order perturbation theory in In terms 

of 



we find 



(2£+l)/2 {2^+l)/2 



Wk{r)dr, 



exp(-i(;i/2) l + 0{t 
l + fe^5 +©(£-). 



(3.17a) 

(3.17b) 
(3.17c) 



We have checked that our numerical solution reproduces the above behaviour for large 1. The expres- 
sion (3.15) for the nucleation rate is finite, as can be easily checked by considering the identity 




1=2 



sinh Dtt 
Z)7r(l +L>2)- 



(3.18) 



This last expression is very useful for the numerical computation. To obtain an accurate result we 



need to compute the exact value of q only up to £ =10-100. Large 
values of k and/or large couplings in the potential. 



> 30 are necessary only for large 



4 Results 

The various steps in our calculation are summarized in fig. 1, for a theory described by the potential 
of eqs. (g), ( p:23cl ) with = -5 • 10"^ fcg, = -6.34-10-2 A:o''^ A^q = 0.1 fco- The dimensionful 
scale /cq is related to the temperature T of the system. It is determined by our assumption that the 
effective three-dimensional description becomes valid at /cq. Previous studies of dimensional reduction 



in the context of the effective average action ]17|, ^ indicate that ~ T. It must be pointed out 
at this point that the two-minimum structure is often a consequence of the integration of fluctuations 
of effective three-dimensional degrees of freedom (additional scalar or gauge fields). In such cases the 
scale ko must be taken sufficiently small for this structure to emerge [27|- p9[| p|. 

In fig. la we present the evolution of the potential Uk{4>) as the scale k is lowered. We have shifted 
the metastable vacuum to (p = 0. The solid line corresponds to k/kQ = 0.513, while the line with 

^The calculation of the pre-exponential factor must take into account fluctuations associated with these additional 
degrees of freedom jjcj] . 
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Figure 2: The behaviour of the nucleation rate I for several values of the parameters of the model with 
^ko/ko = 0.1. We show the values of (diamonds), \n{Ak/k'j) (stars) and — ln{I/kj) (squares) as a 

function of k / J Uj! {(I)t) . 
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longest dashes (that has the smallest barrier height) corresponds to kf/ko = 0.223. At the scale kf 
the negative curvature at the top of the barrier is slightly larger than —k'j. This means that the pole 
of the function L'i{w) at if = —1 is approached in eq. ( 2.1^ ). This is the point in the evolution of 



the potential where configurations that interpolate between the minima start becoming relevant in 



the functional integral that defines the coarse-grained potential |32, For this reason, we stop the 
evolution at this point. The potential and the field have been normalized with respect to kf, so that 
they are of order 1. We observe that, as k is lowered, the absolute minimum of the potential settles at 
a non-zero value of (j), while a significant barrier separates it from the metastable minimum at (j) = 0. 
The profile of the critical bubble (^b(^) is plotted in fig. lb in units oi kf for the same sequence of 
scales. For k ^ kf the characteristic length scale of the bubble profile and 1/k are comparable. This is 
expected, because the form of the profile is determined by the barrier of the potential, whose curvature 
is ~ —k'^ at this point. This is an additional indication that we should not proceed to coarse-graining 
scales below kf. We observe a significant variation of the value of the field (j) in the interior of the 
bubble for different k. This is reflected in the form of the quantity Wk(r), defined in eq. ( |3l9c| ), which 
we plot in fig. Ic. 

Our results for the nucleation rate are presented in fig. Id. The horizontal axis corresponds to 



k/ ^U'l{4)t), i.e. the ratio of the scale k to the square root of the positive curvature of the potential 
at the absolute minimum. The latter quantity gives the mass of the field at the absolute minimum. 
Typically, when k crosses below this mass (corresponding to the value 1 on the horizontal axis) 
the massive fluctuations of the field start decoupling and the evolution of the convex parts of the 
potential slows down and eventually stops. The dark diamonds give the values of the action of the 
critical bubble at the scale k. We observe a strong k dependence of this quantity, which is expected 
from the behaviour in fig.s la-lc. The stars in fig. Id indicate the values of ln(Afc/A;|). Again a 
strong k dependence is observed. More specifically, the value of Ak decreases for decreasing k. This 
is expected, because k acts as the effective ultraviolet cutoff in the calculation of the fluctuation 
determinants in Aj.. For smaller /c, fewer fluctuations with wavelengths above an increasing length 
scale ~ 1/A; contribute explicitly to the fluctuation determinants. The dark squares give our results 
for — ln(//A;|) = — ln{Ak/kf). It is remarkable that the k dependence of this quantity disappears 



as k crosses below yU'i!{((>t) and approaches kf. The small residual dependence on k can be used to 
estimate the contribution of the next order in the expansion around the saddle point. It is reassuring 
that this contribution is expected to be smaller than ln{Ak/kf ). 

This behaviour confirms our expectation that the nucleation rate should be independent of the scale 
k that we introduced as a calculational tool. It also demonstrates that all the configurations plotted 
in fig. lb give equivalent descriptions of the system, at least for the lower values of k. The implication 
is that the critical bubble should not be identified just with the saddle point of the semiclassical 
approximation, whose action is scale dependent. It is the combination of the saddle point and its 
possible deformations in the thermal bath that has physical meaning. 

In fig. 2 we present the calculation of the nucleation rate for several values of the parameters of 
the model. Each row is computed for initial potentials U^q with the same values of fif. and A^q. It 
also has the same value of kf. Thus, it corresponds to the same Z2-symmetric theory in the limit 
7fco = 0- Different values of 7fcp are used for the three calculations in each row, so that the saddle 
points have different profiles and, therefore, different nucleation rates are predicted. Moving down the 
sequence of rows, the value of A^p is kept fixed, while is reduced. The effective dimensionless 

coupling Afco/y^2j^u|j of the Z2-symmetric theory increases, and the resulting potentials have more 
pronounced barriers relative to the location of the minima. This indicates that the effect of fluctuations 
should be enhanced. The last row corresponds to a Z2-symmetric theory that starts approaching the 
Wilson-Fisher fixed point during the evolution of the potential, before deviating towards the phase 



with symmetry breaking or the symmetric one pq , 21, 27, 35]. We have tried to keep the values of 



the predicted nucleation rates comparable in each column of fig. 2. However, the discontinuity in the 
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field expectation value during the first-order phase transition A(j) = (j)t decreases as we move down 
each column. Thus, the strength of the phase transition diminishes. 

The most striking aspect of the comparison of the results in each column concerns the relative 
values of Sk and ln(Afc/A;j). In the first column the contribution of the prefactor to the nucleation 
rate is much smaller than that of the action of the saddle point. The main role of the prefactor 
is to remove the /c-dependence from I/kj- As we move down each column, the difference between 
Sk and ln{Ai^/kj) diminishes. In the last row the two quanitities are comparable. This confirms 
our expectation that the effects of fiuctuations should be enhanced in more weakly first-order phase 
transitions. The second observation concerns the ^-dependence of the predicted nucleation rate. In 
the last row the contribution from the prefactor fails to cancel completely the k dependence of the 
action of the saddle point. In more quantitative terms, when the prefactor is taken into account the k 
dependence of the nucleation rate is reduced by a factor ~ 10 in the first row, while only by a factor 
~ 2 in the last one. The reason for the above behaviour is clear. In the last row and for k = kf, 
the nucleation rate is roughly equal to or smaller than the contribution from the prefactor. Thus, 
the effect of the next order in the expansion around the saddle point is important and can no longer 
be neglected. This indicates that there is a limit for the validity of Langer's picture of homogeneous 
nucleation . For sufficiently weakly first-order phase transitions the saddle point of the semiclassical 
approximation is overwhelmed by the fluctuations around it. As a result, one can no longer rely on a 
picture based on the semiclassical approximation. 

5 Summary and conclusions 

In this paper we addressed the problem of the calculation of nucleation rates for first-order phase 
transitions in the context of high-temperature field theories. The most commonly employed field- 
theoretical tool in such investigations is the effective potential, usually computed within a perturbative 
scheme. A first-order phase transition is expected if two minima are present in the effective potential, 
separated by a barrier. The nucleation rate is calculated through an expansion around the saddle point 
of the functional integral that interpolates between the two minima. The rate is given by an expression 
that involves an exponential suppression by the action of the saddle point, and a pre-exponential factor 
that includes the fiuctuation determinant around the saddle-point configuration. 

Several major obstacles must be overcome before a consistent description can be obtained. The 
most obvious one is a consequence of the convexity of the effective potential. This precludes the 
discussion of tunnelling, as no barrier exists between the minima. All relevant information is washed 
out by the Maxwell construction. In most studies, a perturbative approximation to the generating 
functional of the IPI Green functions is used instead of the effective potential. Such a quantity 
has non-convex parts, but it also has imaginary parts that are difficult to interpret. For radiatively 
induced first-order phase transitions (a common occurrence in field theory) there are more conceptual 
difficulties. It is not clear which fiuctuations of the system generate the two- minimum structure of the 
potential and which are associated with the pre-exponential factor in the calculation of the nucleation 
rate. It is difficult to resolve in a clear way the issue of double-counting the effect of fluctuations. At 
the technical level, a serious issue concerns the ultraviolet divergences that appear in the calculation 
of the pre-exponential factor. They must be cancelled by counterterms of the original action, in a way 
consistent with the calculation of the potential that determines the action of the saddle point. 

In this work we followed an approach that resolves all the above issues. It relies on the introduction 
of a coarse-graining scale k in the problem. This scale separates the high-frequency fluctuations of 
the system, which may be responsible for the presence of the second minimum through the Coleman- 
Weinberg mechanism, from the low-frequency ones which are relevant for tunnelling. The appropriate 
tool for the calculation of the nucleation rate is the coarse-grained free energy at a non-zero value of 
k. The pre-exponential factor is well-defined and finite, as the scale k acts as an ultraviolet cutoff 
in the calculation of the fiuctuation determinants. This is a natural consequence of the fact that all 
fluctuations with typical momenta above k are already incorporated in the form of the coarse-grained 
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free energy. 

We employed the formalism of the effective average action fl^, 18 1, which we identified with 



the coarse-grained free energy. As a starting point, we considered for a real scalar field at a scale 
feo below the temperature, such that the theory has an effective three-dimensional description. We 



approximated by a standard kinetic term and a potential with two minima given by eq. (2.21). 
We assumed that this form of the potential results from the bare potential Ua after the integration 
of (quantum and thermal) fluctuations between the scales A and fco- Some of these fluctuations may 
correspond to additional massive degrees of freedom that decoupled above the scale kQ. If the two- 
minimum structure is a consequence of the integration of fluctuations of effective three-dimensional 
degrees of freedom (additional scalar or gauge fields) the scale kQ must be taken sufficiently small for 
this structure to emerge ||27|]- [^ ]. 

We computed the form of the potential Uk at scales k < k^ hy integrating an evolution equation 
derived from an exact flow equation for F^. Uk is non-convex for non-zero k, and approaches convexity 
only in the limit /c — > 0. The nucleation rate must be computed for k larger than the scale kj at which 
the functional integral in the definition of Uk starts receiving contributions from field configurations 
that interpolate between the two minima. This happens when —k'^ becomes approximately equal to 
the negative curvature at the top of the barrier js^]. For k > kf the typical length scale of a thick- wall 
critical bubble is ^ 1/k. 

We performed the calculation of the nucleation rate for a range of scales above and near kj. We 
found that the saddle-point configuration has an action 5^ with a significant k dependence. For 
strongly first-order phase transitions, the nucleation rate / = Akexp{—Sk) is dominated by the expo- 
nential suppression. The main role of the prefactor A^, which is also /c-dependent, is to remove the 
scale dependence from the total nucleation rate. The implication of our results is that the critical 
bubble should not be identified just with the saddle point of the semiclassical approximation. It is the 
combination of the saddle point and its possible deformations in the thermal bath (accounted for by 
the fluctuation determinant in the prefactor) that has physical meaning. 

For progressively more weakly first-order phase transitions, the difference between and In^A^/kj) 
diminishes. This indicates that the effects of fluctuations become more and more enhanced. At the 
same time a significant fc-dependence of the predicted nucleation rate develops. The reason for the 
above deficiency is clear. When the nucleation rate is roughly equal to or smaller than the contribution 
from the prefactor, the effect of the next order in the expansion around the saddle point is important 
and can no longer be neglected. This indicates that there is a limit for the validity of hanger's picture 
of homogeneous nucleation Q. For sufficiently weakly first-order phase transitions the saddle point 
of the semiclassical approximation is dominated by the fluctuations around it. Despite the presence 
of a discontinuity in the order parameter, one can no longer rely on a picture based on the semi- 
classical approximation. An alternative picture must be developed for the description of the physical 



system |l41 |. 



Finally, we point out that our results are relevant for the question of quantum tunnelling in the 
zero-temperature three-dimensional theory. Only a small numerical factor differentiates between our 
expressions and the ones that determine the quantum-tunnelling rate in that case [^, ^ . The essential 
qualitative conclusions remain unaffected. 
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